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Abstract 

We generalize the relativistic flux tube model for arbitrary two- or three-body systems. The 
spin-independent and spin-dependent contributions of the flux tube to the total Hamiltonian are 

o ' 

£SJ ■ computed in perturbation. In particular, we show that the spin-dependent part exhibits a universal 

spin-orbit form: It does not depend on the nature of the confined particles. The general equations 
we present, being well-defined for light particles, can thus be applied to usual as well as exotic 
hadrons such as hybrid mesons and glueballs. 
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I. INTRODUCTION 



A successful way of understanding the mesons is to approximate the gluon exchanges 
between the quark and the antiquark by a straight string, which is responsible for the 
confining interaction. The relativistic flux tube model (RFTM) is an effective QCD model 
based on this picture Apart from mesons, it has been generalized to baryons and to 
more exotic particles like glueballs and gluelumps (gluon attached to a point-like qq pair) 
j^]. As the present application domain of the RFTM has exceeded its original formulation, 
it is interesting to explicitly write its equations for general two- or three-body systems, in 
order to apply it to some cases of current interest like hybrid mesons seen as qqg states, 
three-gluon glueballs, etc. It is done in Sec. HIl where it is assumed that the dynamical 
contribution of the flux tube is small enough to be treated in perturbation. This approach, 
that we previously called the perturbative flux tube model (PFTM), is rather satisfactory 
since it reproduces the exact RFTM spectrum up to 5% J, |5|]. 

The RFTM and the PFTM neglect the spin of the quarks. Recently, an attempt to 
exactly include spinning particles in the RFTM has been made jfj, but it quickly leads to 
very complicated equations. In this work, we rather compute these spin contributions within 
the PFTM framework in order to get more tractable expressions. Similar calculations have 
already been performed for mesons and baryons with the Wilson loop technique and the 
background perturbation theory [3, Q]. But, the originality of our approach, inspired by 
Ref. [9(, is that the spin correction we obtain remains valid for particles of arbitrary spin 
(not only 1/2 as for the quarks) as well as of arbitrary mass (not only for heavy quarks, but 
also for light and massless bodies). It has a universal spin-orbit form, which is computed in 
Sec. IIHI Some conclusions are finally drawn in Sec. IIVI 



II. PERTURBATIVE FLUX TUBE 



A. Two-body systems 

When the orbital angular momentum i is equal to zero, the two-body RFTM reduces 
to the spinless Salpeter Hamiltonian (SSH) Hq = a/ pf + m\ + m\ + ar, where the 

linear confining potential ar is generated by the straight flux tube linking the quark to the 
antiquark a being the energy density of the flux tube. For not too large £, i.e. i < 6, 
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the dynamical contribution of the flux tube can accurately be treated as a perturbation of 



Hq. The auxiliary field technique allows to compute this perturbation, which is given by 

fly, HQ, 

AH = al(l + 1) [MjJ% + (4 ~ ^1^2) + Qi + fJ-2)ar] , . 

2/ii/i 2 r [12/ii// 2 + 4(/ii + /i 2 )ar + (ar) 2 ] 
In this equation, //j can be interpreted as the constituent mass of the particle whose bare 

mass is rrii. It reads 

' N (2) 



Pi + rn$ 



the average being computed with the eigenstates of the unperturbed Hamiltonian H . Let 
us notice that /ij > even if m ; = 0, as for u, d quarks and gluons. It ensures that the 
correction ([T]) will be defined in every case. In the heavy quark limit where /ij ~ m 8 ^> ar, 
the flux tube contribution becomes 
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in agreement with the spin-independent correction to potential ar as predicted by the Wilson 
loop technique Our formula ([I]) actually generalizes Eq. (jHJ) to the case of light quarks. 
It is worth mentioning that, although the discussion we made was based on the mesonic 
case only, the PFTM Hamiltonian, give n by H = H + AHf t , can be applied to successfully 
describe glueballs and gluelumps [3l. 111]. 



B. Three-body systems 

Lattice QCD calculations support the idea that the so-called Y-junction is the more 
realistic configuration of the static color field in baryons [3]. In this scheme, each quark 
generates a flux tube, the three flux tubes meeting at the point Y which minimizes the total 
energy a^. a« and are respectively the energy density and the length of the flux tube 
starting from particle i. 

We assume here that the value of aj is given by the Casimir scaling hypothesis [l3|] . Under 
that assumption, the energy density aj of a flux tube is proportional to its quadratic SU (3) 
Casimir operator. Then, one should have a g = (9/4) a q . If the three particles are of the same 
nature (quark, gluon, . . . ), the aj are all equal, and Y is the Toricelli point, minimizing the 
total length of the three flux tubes. However, the three bodies can be different in general 
systems. For example, in a hybrid meson seen as a qqg bound state, one can show that, 
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assuming the Casimir scaling hypothesis, the Y-junction is fixed on the gluon P, [ijj], with 
two flux tubes linking the gluon to the quark and to the antiquark. 

The equations defining the three-body RFTM in the center of mass (CM) frame read 

tt 
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(4b) 



The different symbols appearing in these relations have to be clarified. Let xy and vy be the 
position and the velocity of the Y-junction respectively; similarly, X{ and v\ are the position 
and the velocity of particle i. Then, fi = af, — xy, and a« is the energy density of the flux 
tube linking Y to that particle. Moreover, v\ = 9vi± + (1 — 9)vyi±, where Vi± and vyn_ 
are the components of tfj and Vy orthogonal to r*j. In order to find the PFTM Hamiltonian 
corresponding to Eqs. (jl]), one should try to apply the auxiliary field formalism as in the 
two-body case. However, this procedure is too complex here because we are dealing with 
a three-body problem. What can be done is to neglect the string contribution in Eq. (14aj) . 
This approximation leads to Vi = Pi/fii, where ^ is again defined by Eq. (T5]). Then, a 
development of Hamiltonian (14b j) at the order (vf) 2 leads to H = H + AH{ t , where H is a 
three-body SSH with an Y-junction potential, 



ff» = E 



and where 



AH {t = 
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i=i 
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Pi ± -.2 , Pi± VYi± 



(5) 



(6) 
reduces 
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is the dynamical contribution of the flux tubes. In the limit of heavy quarks, Eq. 
to the results of Ref . . 

This term can be further simplified in two cases, following the position of the Y-junction. 
If xy Xi, as in baryons and three-gluons glueballs, we can assume in good approximation 
that the Y-junction is located at the CM. This approximation only overestimates the po- 
tential energy of the genuine junction by about 5% in most cases 15j]. Then, vy ~ as we 
work in the CM frame. By definition, \pi±\ = |-£i|/Vj with Lj = fj x p^ the orbital angular 
momentum for the particle i, and Eq. becomes 



where Li and are now relative to the CM. It is worth mentioning that this last formula 
is a three-body generalization of Eq. ([3]): If the sum in Eq. ([7]) is performed for a two-body 
system only, formula is recovered. 

If the Y-junction is located on one of the three bodies, as it is the case in a hybrid meson 
[sj], we can arbitrarily set xy = £3. Then, r 3 = and vy = V3. Equation ([6]) can now be 
written as 

2 r / -»2 -»2 -> -> \ 1 

(8) 
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This expression is more complicated than (ED) since it involves P3i±, that is the component 
of P3 which is orthogonal to = x.- t — x%. Let us notice that the general equation (jSJ) agrees 
with the results of Ref. jjjj] in the limit of static quarks. 

Let us note that, when the quarks are not static, it has been shown in Refs. [16fl that 
the Y-junction becomes unstable at the classical level. Two quarks actually tend to form a 
diquark linked to the third quark by a single flux tube. Whether such an instability remains 
or not at the quantum level is still a matter of research. Moreover, the relevance of such an 

n 

instability for resonances is questionable [16|. It is also worth mentioning that the Y-picture 
can give good description of baryon spectra, even for highly excited states [l^] • We suggest 
that the predictions of the Y-junction and of the diquark-quark pictures could be compared 
in a quantized version of our perturbative flux tube model. We leave such a comparison for 
future works. 



III. SPIN-DEPENDENT FLUX TUBE CONTRIBUTION 

Let us consider a point-like particle of mass m and charge q, evolving in the vector 
potential generated by a fixed source. Its equations of motion are different following its 
spin, but, in every case, the momentum p^ has to be replaced by 7r M = p^ — qA^ in order 
to take into account the interaction with the external field. In this section, we will focus 
on particles whose spin is either 1/2 or 1 (quarks or gluons). At the quantum level, the 
corresponding equations can all be written in a Schrodinger-like form id t ip = 'Hip-, where ip 
is a "spinor" , whose number of components is different following the spin of the considered 
particle: 4 for a Dirac spinor, and 6 for a spin-1 particle. In this last case, one obtains the 
so-called Duffin-Kemmer-Petiau equation 181 ] . 
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Once Hamiltonian 7i is known, its nonrelativistic limit can be computed thanks to a 



Foldy-Wouthuysen transformation 



191 ] . For a spin-1/2 particle, one obtains for the positive- 



energy part of the Hamiltonian 191 ] 

H = m - qA + ^- - ^-SB + -^(vf x E) + ^V£. (9) 
2m m 2rrr errr 

The symmetrizations of non commuting operators were not written in order to simplify the 
notations. S = a/2 are the spin-1/2 matrices, and the electric and magnetic fields are given 
by E = — VAo, B = Vxi. The Foldy-Wouthuysen Hamiltonian for a particle of spin 1 
has been computed in Ref. [2(|. Its positive-energy part is 

H=m - qAo+ *L-jLaa, do) 

with (S^) = —ie^ the spin-1 matrices. 

Hamiltonians © and ([TO]) were computed in the reference frame of the static source 

under the implicit assumption that was an abelian vector potential. So, they could be 

seen as part of an effective model of QED, but not of QCD. However, this approach can be 

applied to QCD by an appropriate choice of A^ as it has already been shown in Ref. jsj]. We 

will recall here the main points of this work in order to be self-contained. The question is: 

Can the PFTM be simulated by a particular form of the vector potential? The straight flux 

tube of the PFTM actually mimics the configuration of the chromoelectric field as it can be 

□ 

observed in lattice QCD calculations [12J. Consequently, in the rest frame of the flux tube, 
denoted hereafter as the FT frame, the corresponding gluon field is purely chromoelectric. 
The Faraday tensor in the FT frame is thus of the form 



F F J = (W*r - Mr) E a (r) ^, (11) 



with r the distance between the static source and the test particle. The tensor ( TTTT) is such 
that 

K T , FS] = *F bc F a F J u F b Z y « f abC FaE b A c = 0. (12) 

Equation (1121) actually shows that the QCD Faraday tensor behaves as a QED one in this 
special case. Consequently, we just have to search for the most general vector potential 
ensuring E FT = —af/r, B FT = 0, which corresponds to the linearlv rising potential A FT = 
ar giving the static energy of the flux tube. It is shown in Ref. [9J] that this condition is 
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fulfilled if, in the static source's frame, which is also the CM frame, 



ar\ 1 + 
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where p and r are the momentum of the particle and the flux tube length in the CM frame. 
In the nonrelativistic limit, ar Cm, \p \/m <C 1, and thus |9|] 



A t ^ - 

A ps ar, A ps — p, 
m 



^ r -> aL 
E=-a~, B = — . 
r mr 



(14) 



Formula (|14|) can be injected in the Hamiltonians @ and ffTUl) . We set g = 1 because it 
can be absorbed in a redefinition of a. Finally, the spin-orbit terms we were looking for can 
be collected. For a spin- 1/2 particle, 
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in agreement with Ref. 
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For a spin-1 particle 



S(n x E) 



aSL 
2m 2 r ' 
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We can conclude from this discussion that if one takes into account the spin of the 
particles, a spin-orbit correction must be added to the PFTM. It has a universal spin-orbit 
form for spin-1/2 and spin-1 particles, as it can be seen by inspection of Eqs. (TT5"]) and (I16j) . 
and can be thought as a Thomas precession term in the color magnetic field. This relativistic 
correction is consequently given in the case of a general two-body system by 
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in agreement with the Wilson loop formalism in the heavy quark case Q. The last term 
of this equation is only valid in the CM frame, with r the total flux tube length and L the 
relative orbital angular momentum. Formula (|17[) is clearly not valid for light particles (u, d 
quarks, and gluons). However, another approach leads to corrections which are defined even 
in the massless case, that is the background perturbation theory |8j. Within this formalism, 
it is shown that the spin-dependent corrections in mesons are given in the CM frame by 



AH. 



SD 



a 2dV x (r) 
r r dr 



LS\ LS2 
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(18) 
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Vi(r) is a complicated function of the quarks correlators which, at large r, becomes 
Vi(r)\ r _„ = -ar. Then, 



SD 



a 
2r 



Mi 



^2 



(19) 



Equation (|T9|) is equal to (|T5l) up to a substitution of by the dynamical quark masses 
f^i ©, which are always nonzero. In agreement with our results, the correction (fT9"j) also 
holds for gluons jg]. 

For a system in which the flux tubes meet at the CM, like a meson, a baryon or a glueball, 
the spin-orbit contribution is then be given by 



AH so = 



(LiLiSi 



(20) 



with r, the distance between the junction point, located at the CM, and the particle i. 
Li is the orbital angular momentum of particle i, relative to the CM. It is worth noting 
that our formula (I20p agrees with the results of Wilson loop technique jj] and background 
perturbation theory [8j although it has been established in a totally different way. This 
achieves to give us confidence in our result. Let us also remark that the spin-orbit term (|2*U|) 
is expected to be particularly relevant for the glueballs, since the spin interactions are very 
important in this case, but also because of the Casimir scaling of the energy density of the 
flux tubes. 

Formula (I20p has to be modified when the junction is located on one of the particles, say 
particle 3, as for hybrid mesons. The spin-orbit correction is then given by 



AH** = £ 



i=l,2 



O'i £>i Si 



(21) 



with fi = Xi — X3 and £j = fi x pi, in agreement with Refs. [3] for heavy particles. 

Other studies reveal that the Foldy-Wouthuysen Hamiltonian for a particle of arbitrary 
half-integer spin has the same formal structure that the one for a spin- 1/2 particle 21|, and 



that the same conclusion holds for arbitrary integral spin particles [22j. So, we can assume 
that the spin-orbit term we found is formally valid for any spin, even if the spin-1/2 and 1 
are the most relevant cases for our purpose. 



S 



IV. CONCLUSION 



We have explicitly written the equations ruling a general two- or three-body relativistic 
flux tube model, where the dynamical contribution of the flux tube is seen as a perturbation. 
The unperturbed Hamiltonian is a spinless Salpeter Hamiltonian with a linear confining 
potential, while the flux tube perturbation is compatible with spin-independent relativistic 
corrections arising from other effective approaches j^, 8]. Moreover, we have computed the 
spin-dependent part of the flux tube contribution. It appears to have a universal spin-orbit 
form, which does not depend on the spin of the confined particles. 

The perturbative flux tube model we presented is expected to reproduce the exact spec- 
trum of the relativistic flux tube model with an accuracy of about 5% 0, 5]. We think 
that this approach, supplemented by appropriate short-range potentials, is an interesting 
framework to build effective models describing usual [17J as well as exotic hadrons. Indeed, 
our relativistic corrections are more general than the usual ones because they are also valid 
for light particles. We leave the computation of hadron mass spectra with such a model for 
future works. 
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